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Abstract 

One of the most troublesome contributions to the NNLO QCD corrections to 5 — > Xg'^ 
originates from three-loop matrix elements of four-quark operators. A part of this contribution 
that is proportional to the QCD beta-function coefficient (3q was found in 2003 as an expansion 
in rric/mh. In the present paper, we evaluate the asymptotic behaviour of the complete contri- 
bution for rric ^ rnf,/2. The asymptotic form of the /9o-part matches the small-m^ expansion 
very well at the threshold = mh/2. For the remaining part, we perform an interpola- 
tion down to the measured value of m^, assuming that the /^o-part is a good approximation 
at rric = 0. Combining our results with other contributions to the NNLO QCD corrections, 
we find B{B A,7) = (3.15 ± 0.23) x 10"^ for > l.Q GeV in the S-meson rest frame. 
The indicated error has been obtained by adding in quadrature the following uncertainties: 
non-perturbative (5%), parametric (3%), higher-order perturbative (3%), and the interpolation 
ambiguity (3%). 



1 Introduction 

The decay B — > Xg'j is a well-known probe of new physics at the electroweak scale. The 
current world average for its branching ratio with a cut > 1.6 GeV in the S-meson rest 
frame reads [1] 



where the first error is combined statistical and systematic. The second one is due to the theory 
input on the shape function. The third one is caused by the b ^ d'y contamination. 

The total error in Eq. (11. ip amounts to around 7.4%, i.e. it is of the same size as the expected 
(9(a^) corrections to the perturbative transition b —>■ Xf^'^^°'^'-f. On the other hand, the relation 



holds up to non-perturbative corrections that turn out to be smaller (see Section [Tj). 
Consequently, evaluating the Next-to-Next-to-Leading Order (NNLO) QCD corrections to 
b — Xf^^^°^'~f is of crucial importance for deriving constraints on new physics from the mea- 
surements of S — >■ Xg'j. 

In the calculation of 6 — XJ"^''*°°7, resummation of large logarithms (a<j In M^/m^)*^ is 
necessary at each order in as, which is most conveniently performed in the framework of an 
effective theory that arises from the Standard Model (SM) after decoupling the heavy elec- 
troweak bosons and the top quark. The explicit form of the relevant effective Lagrangian is 
given in the next section. The Wilson coefficients Ci{fi) play the role of coupling constants at 
the flavour- changing vertices (operators) Qi. 

The perturbative calculations are performed in three steps: 

(i) Matching: Evaluating Ci(/io) at the renormalization scale fio ~ Mwyfnt by requiring 
equality of the SM and effective theory Green's functions at the leading order in 
(external momenta) / {Mw, fnt). 

(ii) Mixing: Calculating the operator mixing under renormalization, deriving the effective 
theory Renormalization Group Equations (RGE) and evolving Ci{fi) from /io down to the 
low-energy scale fib ^ mb- 

(iii) Matrix elements: Evaluating the on-shell b Xf'^'^^°^'-f amplitudes at fib ~ i^b- 

In the NNLO analysis of the considered decay, the four-quark operators Qi, . . . ,Qq and the 
dipole operators Qi and Qs must be matched at the two- and three-loop level, respectively. 
Three-point amplitudes with four-quark vertices need to be renormalized up to the four-loop 
level, while "only" three-loop mixing is necessary in the remaining cases. The matrix elements 
are needed up to two loops for the dipole operators, and up to three loops for the four-quark 
operators. 

The NNLO matching was calculated in Refs. [2,3]. The three-loop renormalization in the 
{Qi, . . . ,Qe} and {Q7,Q8} sectors was found in Refs. [4,5]. The results from Ref. [6] on the 
four- loop mixing of Qi, • • • , into Q-j will be used in our numerical analysis^] 

^ The small effect (—0.35% in the branchmg ratio) of the four-loop mixing [6] of Qi, . . . ,Qe into Qs is 
neglected here. It was not yet known in September 2006 when the current paper was being completed. 




(1.1) 



r(5^x,7) ~ r(6 ^ '^*°"^7) 



(1.2) 
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As far as the matrix elements are concerned, contributions to the decay rate that are propor- 
tional to |C7(/ih)p are completely known at the NNLO thanks to the calculations in Refs. [7,8]. 
These two- loop results have recently been confirmed by an independent group [9,10]. Two- and 
three-loop matrix elements in the so-called large-/9o approximation were found in Ref. [11] as 
expansions in the quark mass ratio rric/mb- Such expansions are adequate when rric < 'mb/2, 
which is satisfied by the measured quark masses. Finding all the remaining ( "beyond-/5o" ) con- 
tributions to the matrix elements is a very difficult task because hundreds of massive three-loop 
on-shell vertex diagrams need to be calculated. 

In the present work, we evaluate the asymptotic form of the mc-dependent NNLO matrix 
elements in the limit rric S> mb/2 using the same decoupling technique as in our three- loop 
Wilson coefficient calculation [3]. We find that the asymptotic form of the /^o-part matches 
the small-rric expansion very well at the cc production threshold = mh/2. The same is 
true for the Next-to-Leading Order (NLO) matrix elements. Motivated by this observation, 
we interpolate the beyond-/3o part to smaller values of rric assuming that the Po-pait is a good 
approximation at rric = 0. Combining our results with other contributions to the NNLO QCD 
corrections, we find an estimate for the branching ratio at 0{al). 

Our paper is organized as follows. In Section [2|, we introduce the effective theory and 
collect the relevant formulae for the B Xg'-f branching ratio. The contributions that are 
known exactly in rric are described in Section [3l Expressions for the NNLO matrix elements 
in the large-/3o approximation and in the rUc ^ rnb/2 limit are presented in Sections H] and [5l 
respectively. Section [6] is devoted to discussing the interpolation in rric. Section [7] contains 
the analysis of uncertainties. We conclude in Section [HI Our numerical input parameters are 
collected in Appendix A. Appendix B contains a discussion of the cc production treatment in 
the interpolation. 



2 The effective theory 

Following Section 3 of Ref. [12], the 13 Xg'j branching ratio can be expressed as follows: 



B[B ^ X,7], 



B[B Xccv] 



cxp 



Vr 



cb 



6Q!p 



TlC 



[P{Eo) + N{Eo)] , 



(2.1) 



where Ofem = «em(0) — 1/137.036 and N{Eo) denotes the non-perturbative correction]^ The 
TTic-dependence of -B — > XcCu is accounted for by 



C 



ub 



cb 



T[B ^ Xccu] 
T[B^X^euy 



with neglected spectator annihilation. P{Eq) is given by the perturbative ratio 



\Vcb/Vub\^ r[b ^ X^eu] 



VZVtb 



Vc 



cb 



TT 



PiEo). 



(2.2) 



(2.3) 



2 See Eqs. (3.10) and (4.7) of Ref. [12]. The corrections found in Eqs. (3.9) and (3.14) of Ref . [13] as well as 
Eq. (28) of Ref. [14] should be included in N{Eo), too. 
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Our goal is to calculate the NNLO QCD corrections to the quantity P{Eo). The denominator 
on the l.h.s. of Eq. (12. 3p is already known at the NNLO level from Refs. [15, 16]. 
The relevant effective Lagrangian reads 



■^efl — '^QCDxQEd('W; d, S, C, b) + 



4Gf 

71" 



i=l 



i=l 



(2.4) 



where 



= 






{sL1^lUL){uL'l^'hL), 


Qi = 


{sLi;r''CL){cLimL), 


Q2 = 




Qs = 




Qa = 


{sLl>J''bL)T.,{qi^T'^<l)^ 


Q5 = 




Qe = 


(sL7m7M2 7M3^"^L) E,(g7^^7'^^7'^^ 


Q7 = 


^^^2mi,{sLa^''bR)F^^, 


Qs = 





(2.5) 



The last term in the square bracket of Eq. ( 12. 4p gives no contribution at the Leading Order 
(LO) and only a small contribution at the NLO (around +1% in the branching ratio — see 
Eq. (3.7) of Ref. [12]). Consequently, we shall neglect its effect on the NNLO QCD correction 
and omit terms proportional to Vut in the analytical formulae below. However, our numerical 
results will include the Vut terms at the NLO. The same refers to the electroweak corrections 
that amount to around —3.7% in P{Eq) [12,17]. 

The quantity P{Eq) depends quadratically on the Wilson coefficient^ 



(2.6) 



1,...,6, 
7, 



where the "effective coefficients" are defined by 

{Ci(/i), for i 

C7(/i) + E?=i%C,(/i), for^ 
C8(/i) + E?=i;2,C,(/i), for^ = 8. 

The numbers i/j and Zj are defined so that the leading-order b ^ s'-f and b ^ sg matrix 
elements of the effective Hamiltonian are proportional to the leading-order terms in G^^ and 
Gg^, respectively [18]. This means, in particular, that Kij = SijSjj + 0{as)- In the MS scheme 
with fully anticommuting 75, y= (0, 0, -|, -f , -f ) and z= (0,0,1,-^,20,-^) [19]. 



(2.7) 



In Eq. (30) of Rcf. [7], K^j was denoted by G,j/G„. 
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In Eq. (12.61) . we have assumed that all the Wilson coefficients are real, as it is the case in 
the SM. Consequently, Kij is a real symmetric matrix. 

Once the MS-renormalizecj^ coefficients Cf^{jj) are perturbatively expanded 



where 



as(/i6) = 



Air 



the expression for P{Eq) can be cast in the following form: 

[P^'\fib) + Pi'\Eo,fib) + P^'\Eo,fib) 



Here, P*^°^ and P^^^^ originate from the tree-level matrix element of Qj 



p^'Kf^b) 

pt\f^b) 



{21 



(2.9) 



(2.10) 



c?'\^.,))+2cr\^^,)cf>\^^,l 

while P^'' depend only on the LO Wilson coefficients C\^^^^ . 
defined by requiring that it is proportional to products of the LO and NLO Wilson coefficients 
only (Cp^c]^)'^'^). 



(2.11) 

(2) 

The NNLO correction P, is 



3 The corrections P^^^ and P^^^ 

The corrections P2^^ and pj^-* are known exactly in m^. In order to describe their content, we 
expand Kij {Eq, fib) in a^ifib) 



The coefficients K^^^ are easily derived from the known NLO results 



K 



(1) 



i7 



7^(1) 



Rer«-^7jr^. + 20«(5), 
182 8 



for i < 6, 



9 9 



-TT 



liT'Lb + 4^\\>{S) 



(1) 

78 



44 



TT ~ ~ o787 Lb + 2< 



^'"^ 9 27 ■ 2 

i^^) = 2{l + 6,,)4\6), 



'78 



for i,j ^ 7, 



(3.1) 

(3.2) 

(3.3) 

(3.4) 
(3.5) 



The evanescent operators are as in Eqs. (23)-(25) of Rcf. [4]. 
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where 



In 



m 



IS 



(3.6) 



The matrix 7^°)'^^ and the quantities r> as functions of 



15 



(3.7) 



can be found respectively in Eqs. (6.3) and (3.1) of Ref. [20]. The bottom mass is renormahzed 
in the IS* scheme [21] throughout the paper. The charm mass MS renormahzation scale fi^ is 
chosen to be independent from fib. For future convenience, we quote r[ 



(1). 



-6 ri'\z) 



1666 , , ^, 80 . 
h 2\a(z) + b(z)] nr. 



(3.^ 



The exact expressions for a{z) and b{z) in terms of Feynman parameter integrals can be found 
in Eqs. (3.3) and (3.4) of Ref. [20]. Their small-m^ expansions up to 0{z'^) read [22,23] 



16 



- - - - 3C(3) 
2 3 ^ 



1 



6 4 



-Li + 



12 



'7 27r2 vr^ 



^+ 4+ 3 



— i^. 
2 



5ij 

4 " 



2;^ + 



'457 _57r2_^^ _5^2' 
2T6^T8'^72 ^~6 ^ 



&(.) 



z 
2 



2 



TT 



7r 



-3 + — - 



27r2 



^3/2 



1 

2 + ' 



2 

2 2L 



0(z^L^),(3.9) 



25 

12 9' 



1376 137, 27r2 

^ + 2L, H 

225 30 ^3 



z + (1 - 2L,);z^ + ( -ij^ + J + z 



0{z'Li) 



(3.10) 



where 



In 2;. 



(3.11) 



The functions (f)\f{6 = 1- 2Eo/ml'^) with i,j e {1,2,7,8} can be found in Appendix E of 
Ref. [12]. The remaining ones (that affect P(1.6GeV) by ~ 0.1% only) can be read out from 
the results of Ref. [24]. In particular. 



«(5) 



^47 



^48 



54 



6 l-6+-dn + 



lim "^^^^ 



12 mc^mt 



27 



(5) 



(5). 



(3.12) 
(3.13) 



Once all the ingredients of K^^^ have been specified, P2^^ and P^'^^ are evaluated by simple 
substitutions to Eq. f l2.6p 

= E Cp'' C'f <^ (3-14) 
Pf = 2 j: Cl'^^'' C^^ (3.15) 

4 The /3o-part of P^^^ 

('2) 

The only NNLO correction to P{Eq) in Eq. fl2.10p that has not yet been given is P2 . For this 
contribution, we shall neglect the tiny LO Wilson coefficients of Q3, . . . , Qq. The NLO matrix 
elements of these operators affect the branching ratio by only around 1% [20]. Thus, neglecting 
the corresponding NNLO ones has practically no influence on the flnal accuracy. 



Let us split Kjj'^ into the /3o-parts K^j^^° and the remaining parts K^j'^'^^'^ 

= A, nj + P., = K^'^ + K\f--, (4.1) 
where Uf stands for the number of massless flavours in the effective theory, and 

Kf'' - -\^A,^ = -\{}^- \-f) ^ + P., (4.2) 

Following Ref. [11], we shall take nj = 5. Effects related to the absence of real cc production 
in 6 — »• XP'^''*°°7 and to non-zero masses in quark loops on gluon propagators are relegated 
to K^j^'^'^^. Thus, the only mc-dependent contributions to K^^^° originate from charm loops 
containing the four-quark vertices Qi and Q2- 

The explicit k\J^^° that we derive from the results of Refs. [7,8,11,15,20] read 

K?r' = Po Re {-^rf (z) + 2 [a{z) + h{z) - - ^^^} + «(5), (4.3) 

i^i?* = -\K?r\ (4.4) 
K2)/.o _ ^ r_3803 _ 46 ,^80^,^,^/8 ,_98X _ 16 



= /^ol"^-^^ +y^(3)+(^-7r^-yjL,--L^| + 40^rW, (4-5) 

= 2(l + <5,,)0lf"W, for z,jVr. (4.7) 
The small-mc expansion of Rer2^'*(z) up to 0{z'^) was calculated by Bieri et al. [11] 

, 67454 1 247r2 4 / „ , 

Rerriz) = (I1280 -I5207r2-I7l7r^- 5760C 3 +6840L^ 

^ ^ ^ 6561 729 1215 ^ sv ; ^ 



6 



UAOtt^L, - 2520C(3)L^ + 120^^ + lOOL^ - SOLy z 

(43 - 12 In 2 - 3L^) z^/^ — (ll475 - 3807r2 + 967r^ + 7200C(3) 



243 ' ' 1215 

- 1 1 lOL, - 15607r% + 1440C(3)L^ + 990^^ + 260L^ - GOlfj z^ 

+ 2^/2 (62471 - 24247r2 - 33264C(3) - 19494L^ - bOAir^L^ 

243 2187 V sv ; ^ ^ 

. , 2464 „ / 15103841 7912 . 2368,, , 

- 5184L2 + 2160L^ ) z^ tt^z^^^ + + tt^ + ( 3 

6075 V 546750 3645 81 ^ 

147038 , 352 „ 88 „ 512 .A , q/^^.x 

+ L, + TT^L, + Ll Lt z^ + O(z^^^Lt). 4. 

6075 243 243 ' 243 ^ ^' ^ 



/,(2)/3o 

/,(2)/3o 



The function (pjj (5) reads 



0rr(5) = /3o 



57^(5)^6 + 4/ dxF(2."/) 

^0 



(4.9) 



where F^'^''^^^ clS cl function of x = 2E^/mb is given in Eq. (9) of Ref. [8]E1 



The remaining functions </)-^^'^°(5) will be neglected in our numerical analysis. It should 



not cause any significant uncertainty for Eq = 1.6 GeV. For this particular cut, the NLO 
functions (j)\f{S) affect the branching ratio by around —4% only, which is partly due to a 

certain convention in their definitions {4>i]\o) = for {ij) ^ (77), and 077''(1) = 0). An 
analogous convention is used for 0-^'**((5). The known ^77^* (5) affects the branching ratio by 
around —0.4% only. If the NLO pattern is repeated, an effect of similar magnitude is expected 
from the other 0-^ °((5). 

As in Eq. gl]), we can split Pf ^ = Pf + Pf^^"""^ and express P^^^* in terms of Kf^^\ by 
analogy to Eq. fIXTD 

p(2)/3o^ ^ ^(O)cff ^(O)cfr ^(2)/3o_ 
jj=l,2,7,8 

The "~" sign is used above only because we skip j = 3, 4, 5, 6 in the sum. 

(2) . 

5 The full correction P2 i^i ^he limit rric ^ m\)l2 

The present section contains the main new result of our paper, namely the asymptotic form 
of P2 in the limit S> 'mh/2. It has been evaluated by performing a formal three- loop 
decoupling of the charm quark in the effective theory, using the method that has been previously 
applied by us to the calculation of the three-loop matching at the electroweak scale [3,25]. Once 
the charm decoupling scale is set equal to yU^, one recovers the asymptotic form of the matrix 
elements in the large limit. Details of this calculation will be presented elsewhere [26]. 



^ The original calculation of ^"(2,"/) g^j^^j several other contributions to the photon spectrum was performed 
in Ref. [27]. 
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Figure 1: Re(a + b) (left plot) and Rer2 (right plot) as functions of nic/mb = \fz. See the text 
for explanation of the curves. 



The mc-dependence of P{Eo) at the NLO is dominated by Re[a(2) + b{z)]. At the two-loop 
level, we find the following asymptotic form of the functions a{z) and b{z) 

T. / N 34 1/5^ 101\ I fir 1393 \ ,^fl\ 

81 81 z V45 ^ 2025/ \189 33075/ \z^J ^ ' 
Ima(z) = l^vr and Im6(2;) = ^in are exactly constant for z > ^. 

(2) 

For the real part of the three- loop function (z) introduced in Eq. (14. 3p . we find 

^ (2), ^ 112, 27650 1 / 38 ,2 572 , 10427 8 2 

Re r\> (z) = -Ll + + + - Ll + tt^ 

^ ^ ' 9 " 243 6561 z V405 ' 18225 30375 135 



1 (^L" - -^^L + ^^^^^^^^ _ +o(- 
z^ V2835 ^ 893025 125023500 405^ / 

The small-mc expansion of this function has been given in Eq. (14. 8p . 

The dotted line in the left plot of Fig. [1] corresponds to the exact expression for Re(a -|- 6). 
The solid line presents its small-mc expansion up to 0{z^). The dashed lines are found from 
the large-TTic expansion including terms up to 0{z""') with n = 0, 1, 2, which is indicated by 
labels at the curves. The solid and dashed lines in the right plot of Fig. [T] present the same 
expansions of Rerg . No exact expression is known in this case. 

The plots in Fig. [1] clearly demonstrate that a combination of the small-mc expansion for 

rric < mi,/ 2 and the large-mc expansion for rric > mfe/2 (even in the n = case) leads to a 

(2) 

reasonable approximation to Re(a -|- b) and Rerj for any nic- Moreover, no large cc threshold 
effects are seen at rric = 'mb/2. These observations motivate us to calculate the n = term in the 
large-mc expansion of pP''''^ and use it in Section [6] to estimate this quantity for rric < iTLb/2. 

The expressions that we have found for the leading terms in the large-mc expansion of K^j^^^™' 
are presented below. The necessary leading terms of K^j^ are easily derived from Eqs. (15.11) 



8 



and f l5.2p . taking into account that only with i,j>2 do not vanish at large z, and that 



are 2;-independent for i,j = 4, 7, 8. 



^{2)rcm 
-"■22 



if, 



(2)rcm 
27 



{2)rem 
28 



^(2)rGm 
-"^17 



j^(2)rcm 



(2)rcm 
77 



^(2)rem 
-'^78 



^(2)rGm 
-"■88 



where 



,2 



218 

243 



208 ^ 1^ 
81 



-"■27 -"77 



27 ^ 

4736 



127 35 



324 27 



1150 



729 



L 



1617980 



-"27 -"78 ^ 



729 
127 



19683 



3 

20060 
243 



35 



— K. 
6 



(2)rcm 
27 



324 27 

5 3 



(1) 



C(3) + -^L, 



0(1 



16 



- -Lb I K. 



'(1) 

^78 



'■88 



'(1) 



1237 
729" 

< 

32 



27 



27 



D 



20 
27^ 



9 

31823 . 428 . 26590,, , 

-TT^ + — vr^ In 2 + — ^C(3) 



224 



Ld- 



628487 628 



160 



+ 



729 
512 

na-r + 

27 



27 



1 

40(f ^"^(5) + O 



729 



405 



-vr 



2720 ^ 256 2r 



50 

y 

50 



:7r^ - -Ld I 



3 
2 



'(1) 

^78 



16 



27' 



112 



364 
243 



+ X 



(2)rem 
78 



o 



-TT 



'(1) 

'■88 



(2)rcm 



88 



o 



K, 



(l)rem 
47 



K. 



(1) 

47 



In 



Ml 



27 



and the "decoupling logarithm" 



Lf) = Lh 



In 



(5.4) 



(5.5) 
(5.6) 



0(-j,(5.7) 



(5.8) 



(5.9) 
(5.10) 
(5.11) 

(5.12) 
(5.13) 



(5.14) 
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The function (p^j^'^'^™ [5) reads 



, {2)rem/ 



l-<5 



.77 \5) = -4 / dx 

JO 

Stt ax 



9 3 



275 



25 In" 5 + ( 4 + 7(5 - 25' +5^ ) In 5 + 7 - ^(5 - 75^ + 45^ - 



(5.15) 



where F^^'"^ F(2:"«) and F^^'"'^) as functions of x = 2E^/mi, are given in Eqs. (7)-(9) of Ref. [8]. 
The terms proportional to ax = «i^'*(/^ = 'Ti^'^) occur in Eqs. (15.91) and (I5.15P because the so- 
called Upsilon expansion prescription [28] is followed here, which means that is expressed 
in terms of ml^ , and ax in the ratio 

1 - — as(/ib)ax + . . . (5.16) 



m|^°'^ 9 

is treated as independent from as(/ife), i.e. it is not included in the order-counting when the 
0{al{lj,b)) terms are neglected. A reader who wants to bypass this prescription and use another 
kinetic scheme for rrib should set ax to zero in Eqs. (15.91) and (15.151) and at the same time replace 
ml^ by ml°^'^. Next, the latter mass should be perturbatively re-expressed in terms of the chosen 
^kmctic gg|- q£ ^YiQ renormalon ambiguity. We have verified that the ~ 0.4% effect of the 
ax-terms on the branching ratio can be reproduced by using = 4.74 GeV at the considered 
order. 

The mc-independent quantities Xj'^^'^^ and Xgg'*'^"™ in Eqs. (I5.10p and (15.1 ip stand for the 
unknown non-/3o contributions from the two-loop matrix element of Qg in the theory with 
decoupled charm (together with the corresponding bremsstrahlung) . They will be set to zero 
in the numerical analysis. Neglecting them can be justified by arguing that the contribution 
of i^Jg^ to P{Eo) is suppressed relative to that of k!^ by \QdCf/Cf\ ~ |. In effect, Kif'^" 
in Eq. (14.61) affects P{Eq) by around 0.1% only. The suppression factor gets squared for Kgg . 
Thus, neglecting Xj'^'^^^ and Xgg'*'^'''^ is not expected to cause any significant uncertainty. 



Using the results of Refs. [7,8] , one easily derives an expression for the rric —>■ limit of K' 



(2)rcm 
77 



that differs from K^^^^^^ by inclusion of the real cc production contributions (see Appendix B) 
W""! _n^ _ (,Ai> . r \ 587708 32651 628 , 



27 81 ^ ' 9 " V 9 9 
+ ^Trax + Wrr' (5) (5.17) 

We are going to use this formula for testing the mc-interpolation prescription in the next section. 
By analogy to Eqs. (13.140 and (I4.10p . we can write p^'^^^^^ as 

p^(2)rcm^ ^ ^{0)cfr ^{0)efr ^(2)rem_ ^^_^g^ 
i,j=l,2,7,8 

The "~" sign is used above only because i,j = 3, 4, 5, 6 are skipped in the sum. 
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6 Interpolation in rric 



In the present section, we are going to estimate p^^^-'''"™ for rric < mh/2 by performing a 
certain interpolation between its large-mc asymptotic form and an assumed value at rric = 0. 
In particular, we shall assume that the large-/3o approximation is accurate at = 0, which 
may be argued for by recalling that no charm mass renormalization effects arise at that point. 
Since such an assumption is obviously a weak point of the calculation, two alternative forms of 
it will be considered: 

(a) PP-"" — > 0, (6.1) 

(b) P^^^ + pP + pP — ^ pP^\ (6.2) 

The difference between these two cases will serve as a basis to estimate the interpolation un- 
certainty. As a cross-check, we shall also consider the case 

(c) pf-- fcf^^)':^^f^'"(z = 0), (6.3) 

which does not rely on the large-/3o approximation but rather on assuming that the "77" term 
dominates at mc = 0. In Appendix B, issues related to the cc production are discussed. 

The functional form of P2^^'^'^^{z) that we are going to use for the interpolation is a linear 
combination 

p^(2)rcm ^ \r^^\z)\'^ + X2Rer^2\z)+X3Reri^\z)+X4Z^Reri^\z)+X5. (6.4) 

If p^^^^''*'™ was dominated by renormalization effects, the last three terms would be sufficient. 
The first two terms are included because otherwise no In^ z would be reproduced at large z. 
The function |r2^''p is non-analytic at rric = and takes into account the pure four- 

quark operator contribution to the squared amplitude. The function 'Rer2'\z) matches the 
remaining In^ z in the large-m^ asymptotics. We use this function because it is the only 
genuine NNLO four-quark operator matrix element that is known for small z. 

The determination of the coefficients Xi is most easily explained in a specific numerical ex- 
ample. Let us choose the renormalization scales as follows: /iq = 2Mvi/, /if, = m\f /2 and /i^ = 
fTici'n^c)- When the central values of the input parameters from Appendix A are used but only 
rric is retained arbitrary, Eq. 05.181) yields 

p^(2)rcm ^ (9.57 + 7.12) 2 + 50.52 In z + 47.44 + C , (6.5) 

where the contribution from Eq. (15. 4p has been singled out in the coefficient at In^ z (the first 
term). This first term is assumed to determine xi. Then the second term in the coefficient 
at \r? z determines x^- Next, xs is determined from the coefficient at Inz. Finally, by 
adjusting x^ and x^, one can simultaneously match the constant term in Eq. (16.51) and satisfy 
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0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 



Figure 2: pj^^''"'^, p^'^^\ p['' and P^'^ as functions of mjmi, = for /io = 2Mw, [ih = 
/2 and jic = rnc{mc). The other input parameters are set to their central values from 
Appendix A. See the text below Eq. (I6.14p . 



one of the requirements ( 16. ip . ( 16. 2p or ( 16. 3p . In our example, the results are 



1.45 



rW(^)r-K'(0)| 



8.01 Re 



15.63 Re 



d 



16.522— Rern '[zj, 
dz ^ ' 



(b) 



1.45 



15.63 Re 



z)\ — r 



(1) 



+ 8.01 Re 

d 



(2)-rr(0) 
,(1)/ 



(2), 



0.89z— Rer^ '(2) +40.15, 
dz 



,(2)rem. 



where 



(2) 



PPiz 



1.45 



.(1), 



15.63 Re 



z]\ — r 



r?)(z)-r^^^(0) 



(1) 



8.01 Re 



r?)(z)-rf (0) 



8.14 z^Rer^^^(z) + 21.53, 



(6.6) 



(6.7) 



(6.8) 



0) ^ —40.15 determines the very last term in the (b) case. 



Substituting the central value of mf.{m^ from Appendix A sets 2 to 20 ~ (0.262)^ ~ 0.0684. 



Then, P. 



(2)/3o 



(2)rem / 
! I 

(2)rcm / 



la) (^0) 



(b) P','^'--{zo) 



P 



(2)rcm 



3.86 and 
-3.00 
10.98 
4.50 

j(2)rem 



Pi"\zo) 



P?\zo) 



0.87, 

14.85, 

8.36. 



(6.9) 
(6.10) 
(6.11) 



The corrections P2 from Eqs. (16.61) . (16.71) and (16. 8p are shown in the left plot of Fig. [2] 



as functions of \/z = rric /rrib. The other NNLO corrections Pf \ Pf and Pf ' are shown in 



,(2) 
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the right plot. For the same parameters, one 

Pf^ ~ -8.45, (6.12) 



pP^'{z) ~ 10.35 Re r^'^(z) -r^'^(O) 



,(2)/ X ^(2), 



9.57 Re 



r^^\z) - r^^\0) +0.71 P^'^^zq) - 3.86, (6.13) 



(2)/3o/ 



PPiz) ~ 17.00 Rea(^) + 16.83 Re6(2) -31.04 pf^^o)- -10.64. (6.14) 

The sohd hnes in Fig. [2] show the small-mc expansions up to 0{z^). The dashed hnes describe 
the leading terms in the large-mc expansions. The dotted lines correspond to exact expressions. 
The (a), (b) and (c) cases of the interpolated p^'^^'^'^"' are indicated in the plot. The two vertical 
dash-dotted lines mark the la range for ■mc{mc)/ml^ . 

As in Fig. [1], the large- and small-mc expansions nicely match at = mi,/ 2. It is only a 
consequence of the properties of rg^^ and Her 2'^ that determine the ^-dependence of all the 
considered functions. 

It should be stressed that the shape of the curves in Fig. [2] is quite sensitive to the choice of 
renormalization scales. For instance, when the charm mass renormalization scale /ic is shifted 
from its default value /ic = mc{mc) to twice this value, the coefficients at z-j^Iier^\z) in 
Eqs. (16. 6p . (16. 7p and (16. 8p change quite dramatically (to 3.30, —12.34 and —5.09, respectively). 
However, the resulting effect on the decay rate is partly compensated by a correlated change 
of rndfic) in the NLO correction P2^\ The numerical relevance of P^^^^° is very much /in- 
dependent, which compensates the /tft-dependence of as{fib) in the NLO correction. The rela- 
tively small value in Eq. (I6.13P indicates that fif, = mP/2 that is used in this section is in the 
vicinity of the so-called BLM scale [29]. The renormalization scale dependence of our results 
will be discussed in more detail at the level of the branching ratio in the next section. 

From the above results for P^ and 

P(°) ^ 0.1396, (6.15) 
Pf^ ~ -1.515, (6.16) 

pP{z) ~ 3.347- 1.826 Re a(z) - 1.652 Re 6(z) ^ Pi^^(^o) ^ 1-151, (6.17) 

that are calculated using the same parameters, we obtain 

(a) P{Eo) ~ 0.1226 + 0(aem,asKfe) ^ 0.1192, (6.18) 

(b) P{Eo) ~ 0.1295 + 0(ae^,a,Kfe) ^ 0.1261, (6.19) 

(c) P(Po) - 0.1263 + 0(aem,«sKfe) ^ 0.1229. (6.20) 

In the second step above, the electroweak and O (Vub) corrections are added according to 
Eq. (3.10) of Ref. [170 and Eq. (3.7) of Ref. [12], respectively. 



^ In P^'^\z) (|6.14p and P2^\z) (|6.17p . we display only the z-dependence that is due to a{z) and b{z). Other 
z-dependent effects that originate from (S) are very small. 

^ We are grateful to the authors of Ref. [17] for providing us with their results in an extended version that 
allows for arbitrarily varying ^0 Suiid fib- 
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Figure 3: Renormalization scale dependence of B{B — > ^'^^7) in units 10"'' at the LO (dotted lines), 
NLO (dashed lines) and NNLO (solid lines). The upper-left, upper-right and lower plots describe 
the dependence on /ic, fib and /iq [GeV], respectively. 



Next, using Eq. 1^ with N{Eo) ~ 0.0031, one finds 

(a) i3[5-.X,7]^^>,,,Gev - 3.02xl0-^ 

(b) X,7]^^>,,6Gcv ^ 3.19xl0-^ 



(c) B[B^XsjU^,, 



.Y>1.6GcV 



3.11 X 10" 



(6.21) 
(6.22) 
(6.23) 



7 Estimating the uncertainties 

The resuhs (a) and (b) given at the end of the previous section are within the range of around 
±3% from their average, and the result (c) is in between them. The spread of the three 
solutions decreases for higher scales fib- Moreover, the left plot in Fig. [2] seems reasonable 
even as a description of an extrapolation rather than an "interpolation with an assumption" . 
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/if, [GeV] 


2 


2.5 


5 


2.5 


^^c [GeV] 


1.5 


mcirric) 


2.5 


(a) 


3.01 


3.06 


3.14 


3.03 


3.10 


(b) 


3.23 


3.24 


3.26 


3.19 


3.33 


(c) 


3.16 


3.15 


3.18 


3.10 


3.21 


[(a) + (b)]/2 


3.12 


3.15 


3.20 


3.11 


3.22 



Table 1: B{B -^s7) x 10'' in various cases for several choices of /i;, and /ic. The matching scale 
is set to 160 GeV, and the remaining parameters are as in Appendix A. 



Therefore, we shall assign a fixed ±3% to the interpolation ambiguityO Unfortunately, such 
a rough way of determining this error may need to remain in the B — > analysis until a 
complete evaluation of the three-loop matrix elements is performed. Calculating them just in 
the rric = case would help a lot. 

As far as the higher-order ((9(a^)) perturbative uncertainties are concerned, renormalizat ion- 
scale dependence is usually used to place lower bounds on their size. The three plots in Fig. [3] 
show the branching ratio dependence on each of the three scales, once the remaining two are 
fixed at the values that were used in the previous section (yU-o = 2Mvy, = /2 and /ic = 
mc{mi^). Dotted, dashed and solid lines describe the LO, NLO and NNLO results. The NNLO 
branching ratio is defined as the average of the (a) and (b) cases. Stabilization of the scale 
dependence with growing order in is clearly seen. It is especially encouraging in the case 
of /ifo because the interpolation assumptions (16.11) and (16. 2p violate the analytic /i{,-dependence 
cancellation even at 0{a1). The /ic-stabilization is more obvious because it is built into our 
interpolation prescription by construction. The perfect /iQ-stabilization arises both due to the 
smaller value of at the matching scale and to the fact that the NNLO matching conditions 
are complete. 

By studying plots like those in Fig. [3] for various choices of the fixed scales, we have convinced 
ourselves that ±3% is a reasonable estimate of the higher order perturbative uncertainty. At 
the same time, we find that 3.15 x 10~^ is a good central value. It is reproduced (as an average 
of the (a) and (b) cases), e.g., for /iq = 160 GeV, /if, = 2.5 GeV and Hc = 1.5 GeV. Tabled] 
contains our results for the branching ratio in various cases and for several particular choices 
of the renormalizat ion scales. 

The parametric uncertainty is the most straightforward one. Our input parameters and their 
uncertainties are listed in Appendix A. The correlation between C and mcimc) is taken into 
account. The remaining errors that are quoted in the right columns of Tables [21 [3] and H] in 
Appendix A are treated as uncorrelated. This way we find the overall parametric uncertainty 
of ±3.0%. 

Finally, the non-perturbative uncertainties need to be considered. It has been customary for 
a long time to quote the known (9(A^/m^) and (9(A^/m^) non-perturbative corrections as the 

^ Various uncertainties will be added in quadrature later. Thus, each of them should be understood as a 
"theoretical la" rather than a strict range. Of course, such statements concerning theoretical uncertainties are 
never well-defined, but one can hardly improve on this point. 
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dominant ones (see Section VI of Ref. [30] for discussion and references). The subdominant 
0{K^/m\) and (9(A^/(m^mb)) ones are known, too [13]. More recently, the 0(a;<(A^/(mb — 
2E{)Y) correction was evaluated [14]. All these corrections are included in N[Eq) in Eq. ( 12. ip 
and cause around 2.4% enhancement of the branching ratio. However, the non-perturbative 
corrections that arise in the matrix elements of Qi^2 in the presence of one gluon that is not 
soft remain unknown. They scale like agA/mh in the limit mc <^ mb/2 and like a^A^/m^ in 
the limit ^ ?Ti;,/2l^ Since < mf,/2 in reality, we consider agA/mh as the quantity that 
sets the size of such effects. In consequence, we assign a ±5% non-perturbative uncertainty to 
our result. This is the dominant uncertainty at the moment. The very recent estimates [31] of 
similar corrections to the Qj-Qs interference term are neglected here. They are smaller than the 
overall uncertainty of ±5% that both the authors of Ref. [31] and us assign to all the unknown 
0{asA/mi,) effects. 



8 Conclusions 

Our final NNLO result 

B[B ^ X,7]^^>,eGev ^ (3.15 ± 0.23) x 10"^ (8.1) 

is obtained for the input parameters listed in Appendix A and for the renormalization scales 
fio = 160 GeV, fib = 2.5 GeV and /ic = 1.5 GeV. The total uncertainty is found by 
combining in quadrature the ones discussed in Section [71 As it is clearly seen in Fig. [3], the 
NNLO corrections significantly improve the stability of the prediction with respect to the 
renormalization scale variation and, in consequence, reduce the total error. Comparing the 
NLO and NNLO /ic-dependence in Fig. [3l one realizes why all the previously published NLO 
predictions had significantly higher central values than the current NNLO one. 

In order to relate our result with Eq = 1.6 GeV to the measurements with cuts at 1.8 GeV 
(Belle [32]) and 1.9 GeV (BaBar [33]), one needs to compute ratios of the decay rates with 
different cuts (see, e.g., Ref. [34]). This is a non-trivial issue because new perturbative and 
non-perturbative effects become important in the endpoint region. A new calculation of such 
effects has recently been completed [35] but its numerical results were not yet available when 
the average in Eq. (11. ip was being evaluated. 

We have chosen Eq = 1.6 GeV as default here assuming that it is low enough for the cutoff- 
enhanced perturbative corrections [14,35] to become negligible. If this turns out not to be the 
case, one should use our results with lower Eq. For this purpose, we giv^ 

S[5^X,7]^^>,.5Gev ^ (3.18 ± 0.23) X 10-^ (8.2) 
S[fi ^ X,7]^^>,.,oev ^ (3.20 ± 0.23) X 10"^ (8.3) 
B[B Xs^U^^^.,G.v ^ (3.22 ± 0.23) X 10"^ (8.4) 
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We thank M. Bcncke for a clarifying discussion on this point. 

These results do not include contributions from tree-level diagrams with four quark operator insertions 
that are suppressed either by |T4bP or by squares of the small Wilson coefficients C3, . . . , Cg [27]. 
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i3[5^X,7]^.,>i.2Gev ^ (3.24 ± 0.23) X 10-^ (8.5) 
S[S^X,7]^^>,iGev ^ (3.25 ± 0.23) X 10-^ (8.6) 
B[5 ^ X,7]^^>,oGev ^ (3.27 ± 0.23) X 10-^ (8.7) 

Note added 

A numerical analysis of cutoff-related perturbative corrections (see Section [8]) became available 
[36] when the current paper was being refereed. The authors of Ref. [36] use our result given 
in Eq. (18. 7p and combine it with their cutoff-related corrections that turn out to suppress the 
branching ratio at Eq = 1.6 GeV by around 3% with respect to Eq. (18. ip . Such an effect is 
of the same size as our higher-order uncertainty, which means that the unknown 0{al) non- 
logarithmic effects can be as important at Eq = 1.6 GeV as the logarithmic ones calculated 
in Ref. [36] . Therefore, we leave the results of the current paper unaltered, yet not excluding 
the possibility of taking the correction from Ref. [36] into account in future upgrades of the 
phenomenological analysis, once other effects of potentially the same size are collected (see 
Section 1 of Ref . [37]). 
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Appendix A: Numerical inputs 

In this appendix, we collect numerical parameters that are used in Sections [6] and [71 The 
default value of the photon energy cut is Eq = 1.6 GeV. In Section[6], the renormalization scales 
were /iq = 2Miy, fib = /2 and /ic = mc{mc). The final central value of the branching 
ratio in Eq. (18. ip is reproduced e.g. for /xq = 160 GeV, fib = 2.5 GeV and /ic = 1-5 GeV. 

The other parameters are displayed in the tables below. Errors are indicated only if varying 
a given parameter within its la range causes a larger than ±0.1% effect on the branching ratio 
(18. ip . In such cases, the effects in percent are given in the right column of the corresponding 
table. 

Table [2] contains the four quantities that determine the overall normalization factor multi- 
plying P{Eq) in the expression (12. ip for the branching ratio. 
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parameter 


etlect on Eq. f]8.ip 


B{B Xeez>)exp = 0.1061 ± 0.0017 [42] 
C = 0.580 ±0.016 [39,41] 
iVZVtb/Vcbl"^ = 0.9676 ± 0.0033 [43,44] 
acm(O) = 1/137.036 [45] 


±1.6% 
±2.8% 
±0.4% 



Table 2: Parameters that determine the overall normalization factor in Eq. (12. ip . 



parameter 


etlect on Eq. f|8.1|) 


Mz = 91.1876 GeV [45] 




asiMz) = 0.1189 ± 0.0020 [45,46] 


±2.0% 


mf,poie = (171.4 ±2.1) GeV [47] 


±0.5% 


Mw = 80.403 GeV [45] 




ml^ = (4.68 ± 0.03) GeV [39] 


±0.2% 


mc(me) = (1.224 ± 0.017 ± 0.054) GeV [40] 


±2.8% 



Table 3: Experimental inputs that are necessary in the calculation of P{Eq) before including the 
electroweak and 0{Vub) corrections. 



Table [3] contains the experimental inputs that are necessary in the calculation of P{Eq), 
before including the electroweak and 0{Vub) corrections. For as{Mz), we adopt the central 
value from Ref. [46] but conservatively use a twice larger error which then overlaps with the 
one given by the PDG [45]. One should take into account that the phase space factor C in 
Tableland mc{m^ in Table [3] are strongly correlated. For this reason, we take both parameters 
as determined in the same global fit to the semileptonic -B-decay spectra [38-40] instead of 
adopting mc{m^ from Ref. [48]. The normalized correlation coefficient amounts to [41] 

^ {C m,) - (C) (me) r -0.97, method A, 

r = ^ < A.i 

acam, \ -0.92, method B, ^ ^ 

where the meaning of the two methods is explained in Section III of Ref. [40] . Denoting the 
uncertainties in the right columns of Tables [2] and [3] by A^., where x is the relevant variable, 
one finds that the combined uncertainty due to C and rric is given by 

^c,m. = ^Al + Al^ + 2F\Ac^mJ. (A.2) 

Since Aq and A^^ are very close in size (by coincidence) and F is close to —1, the combined 
uncertainty is much smaller than it would be in the absence of the correlation. For our final 
result, we adopt F = —0.92 that gives larger Ac,mc (±1.1% rather than ±0.7%). 

The dependence of mc{mc) on as{Mz) in the semileptonic fit can affect mcimc) by up 
to ~ 10 MeV [41], which would translate to ~ 0.5% in Aq,^. We neglect this effect here. An 
analogous correlation for C is also neglected, as it turns out to be even smaller [41]. 
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parameter 


ettect on Eq. (18. ip 


aera(Mz) = 1/128.940 [49,501 




sin^^iy = 0.2324 [511 




Mhws e [114.4, 1941 GeV (95% C.L) [451 


0.3% 


(115 GeV is used as central) 




{V:,Vuh) / {V:yth) = -0.011 + 0.0182 [43,44] 




Ai = (-0.27 ±0.04) GeV^ [39] 




As ~ i (m|. - ml) ~ 0.12 GeV^ [45] 




pi = (0.038 ± 0.028) GeV^ [39,41] 


±0.4% 


P2 = (0.0045 ±0.035) GeV^ [39,41] 


±0.6% 



Table 4: Remaining parameters that are necessary for the electroweak, 0{Vub) and non-perturbative 
corrections. 



i 






1 


-0.8411 


15.278 


2 


1.0647 


-2.124 


3 


-0.0133 


0.096 


4 


-0.1276 


-0.463 


5 


0.0012 


-0.021 


6 


0.0028 


-0.013 


7 


-0.3736 


2.027 


8 


-0.1729 


-0.614 



Table 5: The LO and NLO Wilson coefficients C^' (fib) at fn, = mf/2 = 2.34 GeV. The 
matching scale /io was set to 2M\y in their evaluation. 



Table H] contains the remaining parameters that are necessary for the electroweak, 0{Vub) 
and non-perturbative corrections. Since we treat the phase space factor C as an independent 
input, Ai is needed only for the small cutoff-related non-perturbative correction [14]. 

The parameters pi and p2 are needed for the O {h? /m\, A^/ [rribml)) contributions to N{Eo) 
that are derived from the formulae of Refs. [13,52]. An important subtlety in this calculation is 
that the O [IS? /m\) corrections to F — > XuCu^ logarithmically diverge when rriu — > 0, so long 
as hard gluon interactions with the spectator quark are neglected. When they are included, 
one encounters new non-perturbative matrix elements whose values are rather uncertain [53]. 
Since these effects cancel in the ratio {P{Eo) ± N{Eq))/C, we follow the approach that has 
been used in the calculation of C in Ref . [39] , namely we neglect the spectator effects and set 
In mf,/m„ to zero. At the same time, we neglect the additional uncertainty in C that is caused 
by this arbitrary procedure (see the comments below Eq. (25) in Ref. [39]). Our final result 
(18.11) is the same as if the considered subtlety did not occur at all. In the future, simultaneous 
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calculations of C and N{Eq) should be performed on the basis of the semileptonic fit results of 
Ref. [34], which would constitute an independent test of the the current calculation. 

Although the Wilson coefficients C^^{fib) are derived quantities, we quote for convenience 
their values at Hb = m]f /2. The LO and NLO ones are collected in Table [5l From among 
the NNLO ones, only C^^^^{^b) is relevant here. We find C'^^^^ {m\^ /2) ~ 16.81 including the 
four-loop mixing (Qi, . . . , Qq) — > Qi but neglecting the four-loop mixing (Qi, . . . , Qq) — >■ Qs- 

The analytical solutions to the NNLO RGE for the Wilson coefficients can be found, e.g., in 
Section 3.3 of Ref. [54]. The resulting explicit expressions for Cf'^^^ [nb] in terms of C\^^^^{ijlq) 
and Tj = as{fio)/(^s{fJ'b) are given in Ref. [6]. As far as as{fi) is concerned, we have used the four- 
loop RGE and applied two independent methods for solving it: a numerical one implemented 
in RunDec [55] and an iterative analytical one that includes (tiny) QED effects, too [54]. Our 
final results are independent of which method is used. However, the intermediate quantities, 
for which more digits are presented, slightly depend on the method. 



Appendix B: The ttIc = case and cc production 

The B Xs'~f branching ratio contains no contribution from cc production because events 
involving charmed hadrons in the final state are not included on the experimental side. Thus, 
the evaluation of 6 X|"^''*°'^7 should be performed accordingly. However, such a definition 
of 6 ^ jj^parton^ may break down for — ^ 0, since logarithmic divergences containing Inmc 

7-)(2)rem 

may arise m 

No such divergences are present in the NLO contributions to P{Eq), and in the other than 
p^(2)rem ]\j]\jLQ Q^gg xhcrc, all the logarithms of get multiplied by positive powers of this 
mass. If p^^^-'''"™ turns out to be logarithmically divergent at rric 0, it should be redefined 
to include the cc production contributions. The assumptions fl6.ip - fl6.3l) would then refer to 
the redefined quantity that must be convergent at rric — * 0. At the same time, the large- 
rric asymptotics would remain unaltered. After performing the interpolation of the redefined 
quantity as in Section [6l one should subtract the cc production effects at the measured value 
of rric. At present, possible effects of such a subtraction are understood to be contained in the 
interpolation ambiguity. 

Actually, Eq. fl5.17p does contain the cc production contributions in the case. 

However, it is used only in the assumption (16.31) that serves just as a cross check, and has no 
influence on our flnal numerical result (18. ip . 
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